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C _^ Abstract. We consider a certain finite group for which Kloosterman sums 

£Nj appear as character values. This leads us to consider a concrete family of 

I commuting hermitian matrices which have Kloosterman sums as eigenvalues. 

'~i These matrices satisfy a number of "magical" combinatorial properties and 

5 they encode various arithmetic properties of Kloosterman sums. These matri- 

^v^j ces can also be regarded as adjacency matrices for multigraphs which display 

Ramanujan-like behavior. 

H 

c~! 1. Introduction 

For a fixed odd prime p, let £ = exp(2m/p) and define the classical Kloosterman 
^ sum K(a,b) := K(a,b,p) by setting 

p-i 

(N K(a,b) = Y,( an+m f 1 - 1 ) 

J> n=l 

where n denotes the inverse of n modulo p. From (1.1 1, it follows that K (a, b) is real 
and that its value depends only upon the residue classes of a and b modulo p. In light 
C*~) of the fact that K (a, b) = K(l, ab) whenever p \ a, we focus our attention mostly 

t^j- on Kloosterman sums of the form K(l,u). Moreover, we adopt the shorthand 

K(u) :— K{1, u) or even K u :— K(l, u) when space is at a premium. 

In the years since they appeared in Kloosterman's paper on quadratic forms 
[12) . these exponential sums and their generalizations have found many diverse 
^ applications. We do not attempt to give a historical account of the subject and 

instead direct the reader to [SJ [7J |H1 E] ■ 

In this note we construct a certain finite group for which Kloosterman sums 
appear as character values (Section |2]). This eventually leads us to consider a 
concrete family of commuting hermitian matrices which have Kloosterman sums as 
eigenvalues (Sectionj3l) . These matrices satisfy a number of "magical" combinatorial 
properties (SectionETand they encode various arithmetic properties of Kloosterman 
sums (Section [5]). Moreover, these matrices can be regarded as the adjacency 
matrices for multigraphs which display Ramanujan-like behavior (Section [6]). 
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2. The group G and its representation theory 
Let p > 3 be an odd prime and define the subgroup 



G 



x y \ I x 1 z 



I ^ V 1 



x e (z/pZ) x , y ,z e z/pZ 



of GL^Z/pZ). Here we identify direct sums of two 2x2 matrices with the corre- 
sponding 4x4 matrices. In the following, we denote matrix groups by bold capital 
letters (e.g., G) and their elements by capital letters (e.g., / denotes the 4x4 
identity matrix in G). Elements of Z/pZ are represented by lower-case letters. 
Letting 



N 



i)®(o 1 



y,zeZ/pZ\, (2.1) 



xG(Z/pZ) x , (2.2) 



x 0\ (x- 1 
l) ® ^ 1 

we find that G = NT and N n T = {/}. Since |T| = p - 1 and |N| = p 2 , we have 

\G\ = (p-l)p 2 . 

The conjugacy classes of G are easily computable and are given in Table [2~Tj 

Since the commutator subgroup [G, G] = N of G must belong to the ker- 
nel of any one-dimensional representation tt : G — > C, it follows that tt(NT) = 
h(N)k(T) = ?r(T) for all N e N and T e T. Thus the one-dimensional repre- 
sentations of G correspond to one-dimensional representations of T = (Z/pZ) x . 
Fix a primitive (p — l)st root of unity £. If T denotes a generator of T, then for 
n = 1, 2, . . . ,p — 1 the formula 

7r(T)=f\ 7r(iV) = l, iVeN 

yields p—1 distinct irreducible representation of G. 

Let us now identify the remaining irreducible representations. As before, we let 
£ = cxp(2ni/p). Fixing a,b € Z/pZ, at least one of which is nonzero, we claim that 
the map tt : G -> End(C[(Z/pZ) x ]) defined by 

X y \ ffi ( X 1 Z 11 = ^az(xh)+by(xh) ^2.3) 



1/ \ 1 



for /i £ (Z/pZ) x is an irreducible representation of G. Verifying that tt is a homo- 
morphism is straightforward, so we only prove irreducibility. 



Suppose that a ^ and note that setting x = 1, y — 0, and z — 1 in (2.3) yields 



for h E (Z/pZ) x . The preceding is just another way of saying that 

t((J ?)©(J })) =diag((C) 1 ,(C Q ) 2 ,...,(C a r 1 ) (2-4) 

with respect to the standard basis {5i,S 2 , ■ ■ ■ , £> P -i} of C[(Z/pZ) x ]. Since n ^ 
it follows that ( a is a primitive pth root of unity and hence the diagonal entries 
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Type 1: 

) — 1 classes 



Ci 



C-2 



V 

2y 



o 1 

1 y- 1 
1 



■■ye (z/ P z)> 



(p — 1 elements) 
(p — 1 elements) 



Type 2: 

2 classes 



Cp + i 



Type 3: C p+2 

1 class 



Type 4: 

p — 2 classes 



-P+3 



i (p-i)iA m A jr 1 



o 1 



:w e (i 



1 2/ 

1 

1 

1 

1 
1 

9 v 

o 1 



g 2 y 
o 1 



1 o 

1 

1 y 

1 

1 o 
o 1 



o 1 



:ye (Z/pZ.y 
:ye (Z/pZy 



:y,z£ 



' 2 (\ :y,zez/ p z\ 



(p — 1 elements) 
(p — 1 elements) 
(p — 1 elements) 
(1 element) 
(p 2 elements) 
(p 2 elements) 



C2p 



g p ~ 2 y 

1 



-(p-2) 



: j/, z £ X/pZ > (p elements) 



TABLE 2.1. The conjugacy classes of G (here g denotes a primitive root 
modulo p). In particular, G has a total of 2p conjugacy classes whence there 
exist precisely 2p distinct irreducible representations of G [2] Theorem 27.22]. 



of (2.4) are distinct. Thus the subspaces of C[(Z/pZ) x ] which are invariant under 
the matrix (2.4) are precisely those of the form span(_ftT) for some K C (Z/pZ) x . 
Suppose that K ^ 0, K ^ (Z/pZ) x , and x £ K. For each k E K, implies 
that 

8( x k-i)k =&xi span(JsT). 

as claimed. The 



xk 




x- x k 




Thus span(A') is not invariant under 7r whence 7r is irreducible 
proof in the case b ^ is similar. 

The choices a = 1, b = and a = 0, & = 1 lead us to two special characters, 
whose values on the various conjugacy classes can be found via a geometric series 
argument. We are interested primarily in tt arising when a ^ and 6^0. Let 
Xj denote the trace of tt corresponding to j = a _1 6. Using the transformation 
rules for Kloosterman sums we find that Xji^k) = K(a,bk) = K(l,jk) — Kjk for 
1 < 3-,k < p— 1. Since the Kloosterman sums K (1), K(2), . . . ,K(p— 1) are distinct 
[H Prop. 1.3], it follows that the characters Xi,X2> ■ • ■ iXp-i are distinct. Now we 
can complete the character table for G (Table 2.2). 
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<0 










i-H O O 



O O 




O O 



I I 



I I 



I I 



s, + 



X 



X 



ii a 

CO 

^ s 

II 3 



CD CD 
O > 



P H 



1 M 

a, c 



a, 
-a 

CD 
X 



o 

bjO 

-3 
o 







is 5P 
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3. The main construction 

3.1. The crucial lemma. From the representation-theoretic information com- 
puted in Section [2j we will construct a family of commuting hermitian matrices 
which encode many fundamental properties of classical Kloosterman sums. Our 
primary tool is the following lemma, which is a modification of 1 1 31 Lem. 4] (al- 
though there the reader is simply referred to [21 Section 33] to compose a proof of 
this lemma on their own). We provide a detailed proof for the sake of completeness. 

Lemma 3.1. Let G be a finite group having conjugacy classes C\,Ci, ■■■ ,C S and ir- 
reducible representations 7Ti,7T2, ■ • • ,7T S with corresponding characters Xi)X2i • • • iXs- 
For 1 < k < s, fix z = z(k) £ C k and let Cij.k denote the number of solutions 
(xi,yj) € Ci x Cj of xy = z and then let Mi — (ci,j,k)j k=i f or •"■ — * — s - 
IfW = {wj,k) S j k = \ denotes the s x s matrix with entries 

,„ _ \Cj\Xk(Cj) . , 

Wj.k - —r. , (3.1) 

dim7Tfc 

and Di — diag(wi,i, Wi^, ■ ■ ■ , w%,s)> then W is invertible and 

MiW = WD, (3.2) 

for i = 1, 2, . . . , s. Moreover, if we let Q = diag(y]Ci|, ylCal) ■ ■ • > \/\C s \), then the 
matrices Ti — Q~~ l MiQ are simultaneously unitarily diagonalizable. To be more 
specific, we have TiU = UDi for i = 1, 2, . . . , s where 

U =^{VWi\x k (C 3 )Y i (3.3) 

is a unitary matrix. 

Proof. For j = 1, 2, . . . , s define 



C, = 



\ ^ 



3 - 

xG Cj 



and observe that 



C,C , = ^2 c h3,kC k (3.4) 
fc=i 

holds for 1 < i,j,k < s. Upon applying \k to Cj we also note that 

Xk (Cj) = \Cj\xkiCj) (3.5) 

since the class function \k assumes the constant value Xfc(Cj) on Cj. 

Since each Cj belongs to the center Z(C[G]) of C[G] (in fact {Ci, C 2 , . . . , C s } is 
a basis for Z(C[G]) by Thm. 27.24] or Thm. 2.4]) and each n k is irreducible, 
it follows that nf.(Cj) is scalar for 1 < j, k < s (this follows from a standard version 
of Schur's Lemma [5] Thm. 29.13]). Thus there exist constants Wjk such that 

■K k {Cj) = Wj, k Id h (3.6) 

for 1 < j, k < s where d k = dim nk and Id k denotes the dk x dk identity matrix. 
Taking the trace of the preceding yields 

Xk(Cj) = d k Wj,k- (3.7) 



Comparing (3.5) and (3.7) we find that 

\ C j\Xk(Cj) = d k Wj. k , 
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which gives us the formula (3.1). Applying 7r r to (4.1) and using (3.6) we obtain 



i,rld r w j,rldr — ^ Cjj^Wkrld 



which clearly implies that 



k=l 



k=l 



Now simply observe that the preceding is the (j, r)th entry of the matrix equation 
( [372] ) . Next we note that W = ^/\Q\QUR where R = diag(dj~ \ d^ 1 , . . . , dj 1 ). In 
particular, it follows that 



QUDiR = QURD, 
= MiQUR 



(R, Di are diagonal) 



(by (3.2)) 



whence QUDi — M^QU since R is invertible. Since Q is invertible this yields TJJ = 
UD l where T { = Q~ 1 M l Q. The fact that \G\~ 1/2 U is unitary (whence W is invert- 
ible) follows from the orthogonality of the irreducible characters Xij Xii ■ • • i Xs- D 



3.2. Main construction. We now apply Lemma [3TT] to the group G constructed 
in Section [2] As we shall see in Section [5] the matrices produced encode many of 
the basic properties of Kloosterman sums. 

Recall that the (j, k) entry (Mj)^ of Mj is defined to be the integer Cij t k de- 



scribed in Lemma 3.1 Since G has four distinct types of conjugacy classes (see 



Table 2.1), we partition each Mi into 16 submatrices. As in Table 2.2 we adopt 
the convention that / = p — 1. It turns out that for 1 < i < f each of the Mi has 
basically the same structure as Mi, so we only display Mi explicitly: 



Mi 



Cl,l,l 


Cl,l,2 ' 


• c i,i,/ 








/ 





■• 


' \ 


Cl,2,l 


Cl,2,2 


■ c l,2,/ 


1 


1 








■• 


■ 


C l,/,1 


c l,/,2 ' 


• c i,/,/ 


1 


1 








•• 


■ 





1 


1 





1 








•■ 


■ 





1 


1 


1 











•• 


■ 


1 




















•• 


■ 




















/ 


•■ 


■ 























/ ■• 


■ 























■• 


■ / / 



(3.8) 



We are interested primarily in studying the entries Cij k for 1 < i,j, k < f and we 
discuss them at length in Section [4] 
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In general, Mi for 2 < i < f differs from Mi only in the upper-left 3x3 blocks. 
For instance the upper-left corner of M 2 looks like 



C 2,l,l 


02,1.2 


c 2,1.3 ' 


• C 2,l,/ 


1 


1 





C2,2,l 


C2,2,2 


C2,2,3 


• c 2,2,/ 








/ 


C2,3,l 


C2,3,2 


C2,3,3 


c 2,3,/ 


1 


1 





C 2,/,l 


C 2,/,2 


C 2,/,3 ' 


' C 2,f,f 


1 


1 





1 





1 


1 





1 





1 





1 


1 


1 











1 


















(3.9) 



We therefore restrict our attention mostly to the case i = 1 since the computations 
for i = 2, 3, . . . , / are almost identical. 

Examining the character table (Table 2.2) of G tells us that 



Di = diag(if H , K 2l ,..., K {p _ 1)tl -1, -1, /,/,..., /), 

p — 1 times 

Q = diag( VT) V7 , ■ • ■ , V7 , 1,P>Pj • • • )P), 

p + 1 times p — 2 times 



(3.10) 



and 





K 2 


■ K f 


-1 


-1 


/ 


f 


f 


/ 


K 2 


K 4 


■ K 2f 


-1 


-1 


/ 


f 


f 


/ 


Kf 


K 2f ■ 


■ K f 2 


-1 


-1 


/ 


f 


f 


/ 


-1 


-1 


-1 


/ 


-1 


/ 


f 


f 


/ 


-1 


-1 


-1 


-1 


/ 


/ 


f 


f 


/ 


1 


1 


1 


1 


1 


1 


1 


i 


1 

















v 1 




P 2 H 2 


p 2 e~ 2 

















P 2 


P 2 e 




p 2^2(p-2) 

















P 2 


P 2 e- 2 


p 2^2(p-2) 


. P V?- 2 > 2 



Since Lemma 3.1 guarantees that the matrices Tj := Q~ 1 M i Q are simultaneously 
unitarily similar to the corresponding ZVs, each of which has only real entries, it 
follows from the Spectral Theorem that each Tj is hermitian. For instance, 



Ti = 



Cl,l,l 


CI, 1,2 


• c l,l,/ 








^7 





• 


■ 


Cl,2,l 


Cl,2,2 


■ c l,2,/ 


1 


1 








■ 


■ 


c l,/.l 


c l,/,2 ■ 


■ c l,/,/ 


1 


1 








■ 


■ 





1 


1 





1 








• 


■ 





1 


1 


1 











■ 


■ 






















■ 


■ 




















/ 


• 


■ 























/ • 


■ 























■ 


■ / 
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Moreover, the unitary matrix U of Lemma |3.1| is given by 



U 



V 





K 2 




-1 


-1 


1 


1 


1 


1 


K 2 


K A 


■■ K 2f 


-1 


-1 


1 


1 


1 




K f 


K 2f 


■ ■ K f2 


-1 


-1 


1 


1 


1 


1 


-1 


-1 


-1 


/ 


-1 


1 


1 


1 


1 


-1 


-1 


-1 


-1 


/ 


1 


1 


1 


1 


•/J 


V7 


■■ VI 


VT 


V? 


I 

v7 


1 

■Jl 


1 

^ '- 


1 

- 






















p 

v7 
p 

•/J 


Pi 
v7 
P« 2 

V? 


s/T 

p£ 4 
vT 


p("- j 
✓7 

















p 

v7 


v7 


p 5 2<P-2) 

■ST 


pe<f- 2 ' 2 
v7 



and it has the property that Till = UDi. In particular, the fcth column of U is an 
eigenvector of Tj corresponding to the fcth diagonal entry of Di. 

In light of the block upper-triangular structure of U and the block of zeros in 
the upper-right of Tj, it follows that the equation TiU = UDi still holds if we 
truncate all matrices involved to their upper left (p + 2) x (p + 2) blocks. We do 
so in order to remove entries that are irrelevant for our purposes and contain no 
useful information about Kloosterman sums. Performing this truncation we now 
consider instead the (p + 2) x [p + 2) matrices 



A = diag(K li ,K 2 i 



>-K"( P -i)i,-l,-l,p- 1) 



(3.11) 



and 



Ti = 



Ci, 1,1 


Ci,l,2 • 


• C i,l,/ 










Ci,2,l 


Ci,2,2 • 


' C *,2,/ 


1 


1 







C *,/,2 ' 


' C *J,f 


1 


1 








1 


1 





1 








1 


1 


1 








VJ 














o / 



(3.12) 



Unfortunately, the new U obtained by truncating the original U is no longer unitary. 
However, this can easily be remedied by normalizing the (p + 2)nd column, leading 
us to redefine U as follows: 



U = 





K 2 ■ 


• K f 


-l 


-l 


77\ 


K 2 


K± ■ 


■ K 2f 


-l 


-l 


77 


Kf 


K 2f ■ 


■ K P 


-l 


-l 


77 


-1 


-1 • 


■ -1 


/ 


-l 


77 


-1 


-1 • 


■ -1 


-l 


/ 


77 


77 


77 • 


■ 77 


77 


77 


i / 



(3.13) 



In summary, the truncated matrices (3.11), (3.12), and (3.13) satisfy T,C7 = UDi 
for 1 < i < p — 1 . 
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Fori = 1,2, 



4. SUBMATRICES OF THE Tj 

,p— 1 we let Bi — {ci,j t kfjk-\ denote the upper left (p— 1) x (p— 1) 



submatrix of Tj (3.12 1. In this section we examine the structure of these matrices. 
Some of these properties will be used in Section [5] to study Kloosterman sums and 
in Section [6] to construct Ramanujan multigraphs. 

4.1. Computing the entries. We claim that the entries of the Bi are given by 

(4.1) 



= i 



P 



where 



i 2 + j 2 + k 2 



2ij -2jk-2ik 



(4.2) 



and (^) denotes the Legendre symbol 



— 1 if a is a quadratic nonresidue modulo p, 
= { if a = (mod p), 

1 if a is a quadratic residue modulo p. 

In particular, it follows that Cj ,• j. € {0, 1, 2} for all 1 < i,j,k < p — 
(4.1) and (4.2), we also have 



1. In light of 
(4.3) 



Ci,j,k — C (r(i),cr(j),cr(k) 

for any permutation a of {i,j, k} and 

c i,j,k = c li,lj,lk (4-4) 

for 1 <i,j,k,l < p — 1 (here the subscripts li, Ij, Ik are considered modulo p). Let 
us now justify the formula (4.1 ) for the entries of Bi. 
According to Lemma ,' 



3.1 



the entries Cjj-^ denote the number of solutions (X, Y) € 
Cj x Cj to the equation = Z for some fixed Z £ C}.. For 1 < i, j, k < p — 1 we 
consider the equation 



1 


ix 







1 








1 x' 1 






1 



\ 



i jy 






1 1 


\ 


1 




1 




1 y- 1 






1 


1 




1 J 




V 





1 / 



(4.5) 



xec, 



YeCj 
i _ 



which instantly reveals that x^ 1 + y 
equation ensures that x, y 7^ 1 so that y 
second equation we obtain the quadratic 



1 and ix 

- x(x — 1) 



3V 



zec k 

- k. Note that the first 



+ (J — k — i)x + k 



1 . Substituting this into the 

0, (4.6) 

which has cither 0, 1, or 2 solutions in Z/pZ. Since k 7^ 0, it also follows that every 
solution x to (4.6) belongs to (Z/pZ) x and hence there is a bijective correspondence 
between solutions (X, Y) € Ci x Cj to (4.5) and solutions x € Z/pZ to (4.6). 
Substituting (2i) _1 [a; — (j — k — i)] for x reveals that (4.6) has the same number of 
solutions as 

x 2 = (j - k -i) 2 - Aik = /3(i, j, k) 



where the function (3(i,j,k) is defined by (4.2). This establishes (4.1). 
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Before proceeding, we should remark that the appearance of the preceding qua- 
dratic is not surprising when one considers the well-known formula 



KM 



p-i 

E 

n=0 



4u 



(4.7) 



which can be found in |3J Lem. 1.1], [T51 eq. (1.6)], or [2TJ eq. (51)]. 



4.2. Rows and columns. Fix l<i,k<p~l and note that as X runs over the 
p — 1 elements of Cj, the variable Y = X~ X Z runs over / distinct elements of G. 
Therefore the sum of the fcth column of Mi must equal /. In light of (3.8), we 
obtain the following formula for the column sums of the Bf. 



(4.8) 




By symmetry, the same formula holds for the row sums of Bi 
For 1 <i,j<p—l fixed we have 

P(i,j,k) = k 2 -2(i+j)k + (i-j) 2 , 



(4.9) 



which we now consider as a quadratic in the variable k. By (4.1) it follows that 
c ij,k = 1 if and only if k G (Z/pZ) x is a root of the preceding quadratic. Since 
i and j are fixed, this holds for at most two values of k. Therefore each row (or 
column) of Bi can contain at most two l's. Let us be more specific. 

• Since p — 2 is odd, it follows from (4.8) that the ith row of Bi contains 
exactly one 1. The remaining p 
which add up to p 

and — r 



2 entries of the ith row are 0's and 2's 
3 by (|4.8|). Thus exactly of these entries are 2's 



of them are 0's. 



• Since p — 3 is even, it follows from (4.8) that for j ^ i the jth row of Bi 
contains either zero or two l's. If the jth row contains zero l's, then 

of its entries must be 2's. If the jth row contains two l's, then of its 
entries must be 2's. 

• Now suppose that j ^ i. We claim that if ij is a quadratic residue modulo 
p, then the jth row of Bi contains exactly two l's and zero l's otherwise. 



The only way to obtain a 1 in the jth column of Bi is for (4.9) to be 
congruent to modulo p for some 1 < k < p — 1. Using the substitution 
k *-> i + j + k (mod p), we see that this occurs if and only if 



k 2 = 4ij (mod p) 



(4.10) 



for some k in {0, 1, 2, . . . , i + j — 1, i + j + 1, . . . ,p — 1}. The forbidden value 
i + j poses no problem since if (i+j) 2 = 4«j (mod p), then p\(i — j) whence 
j = i, contradicting our hypothesis that j ^ i. Thus ( |4.10 1 has a solution 
k 7^ i + j if and only if ij is a quadratic residue modulo p. 

Putting this all together, we obtain Table |4.1| which describes the number of 
elements of each type in a given row/column of Bi. Using this data and the fact 
that there are exactly nonzero quadratic residues of the form ij (j ^ i) and 
nonresidues we can compute the total number of 0, 1, 2's in the matrix Bi (see 
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Row # 


#0's 


#l's 


#2's 


3 = i 


p-i 

2 


1 


p-3 
2 




p-1 
2 


2 


p— 5 
2 




p+1 
2 





p-3 
2 



TABLE 4.1. Number of elements of each type in a given row of Bj. By 
symmetry, the same data applies to the columns of Bi. 



#0's #l's 



#2's 



Total II ip(p-l) \ P -2 I i(p-2)(p-3) 



TABLE 4.2. Total number of elements of each type in the matrix Bi. For 
large p the entries are roughly evenly split between O's and 2's (i.e., approxi- 
mately \p 1 ■ On the other hand, the total number of l's in the matrix is only 
of order p. 



Table 4.2). We can also use this information to compute the sum of the squares of 
the entries of Bi (i.e., the quantity tr B*Bi = tvBf): 

trB? = (p-2)+2(p-2)(p-3) 



2p 2 -9p+ 10. 



(4.11) 



4.3. Magical properties. Along the main diagonal of Bi we have j = k so that 



1 + (- — — ). For j = 1, 2, . . . ,p — 1, this yields the sequence 



4i, 



8i, 



i 2 -Ai(p-l) (mod p). 



(4.12) 



Note that the sequence i 2 — Aij = i(i — Aj) cannot assume the value i 2 since p \ Aj. 
On the other hand, — Aj) assumes every other value in Z/pZ exactly once. Thus 



we conclude that appears in the sequence (4.12) exactly once. Therefore exactly 



one of the diagonal entries of Bi is equal to 1. Since Bi is symmetric, it follows 
that there are an odd number of l's among its entries, in agreement with the data 
in Table HH 



The trace of Bi is easily computed using the above. Since (4.12) assumes every 
value in (Z/pZ) x apart from 1, it follows that there are precisely nonzero 
quadratic residues on the list. Since we already know that a single 1 appears on 
the diagonal of Bi it follows that 

tvB t =p-2. (4.13) 
Next we observe that certain "broken diagonals" of Bi also enjoy curious sum- 



mation properties. Indeed, using (4.3) and (4.4) for j and k fixed we have 

p-i 



p-i 

i=i 



1=1 



p-1 

r=l 




P-2 if j = k, 
3 i£j?k, 



(4-14) 



by (4.8). Define an equivalence relation ~ on pairs (j, k) with 1 < j, k < p — 1 by 
setting (ji, ki) ~ (j 2 , k 2 ) if and only if j\ = Iji (mod p) and ki — lk 2 (mod p) for 
some 1 < I < f. This partitions the indices (j, k) into p — 1 equivalence classes 
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Note: These matrices arc 
best viewed in color (each 
broken diagonal is repre- 
sented using a different 
cofor) . 



of p — 1 elements each, the sum over each equivalence class being given by the 
preceding formula. 

Putting this all together we obtain two different "magic matrices" that are nat- 
urally associated to Bi. First observe from (4.8), (4.131, and (4.14) that if we 
subtract 1 from Ci^j we obtain a new matrix B[ for which each row, column, diag- 
onal, and "broken diagonal" sums to p — 3. For instance, if p = 7 and we subtract 
1 from the (1,1) entry of Bi we obtain the 6x6 "magic" matrix 



/ 1 







1 





o\ 





1 





1 







2 







2 








1 


1 


2 





















2 


2 


V 


2 







2 


/ 



each row, column, and broken diagonal of which sums to 4. 

We can also augment the {p — 1) x (p — 1) matrix Bi with one additional row 
and column from the larger matrix Ti to obtain a px p matrix Ai which also enjoys 
"magic square" properties (i.e., Ai is the upper-left p x p principal submatrix of 
Ti). In particular, each row and column of Ai sums to p — 2 while each diagonal 
and "broken diagonal" sums to p — 3. For p = 11, we obtain the 11 x 11 "magical" 
submatrix A\ of T\ 



( 








1 


2 










2 


2 


\ 





2 


2 


2 





2 














1 





2 


1 





1 


2 




2 








1 


1 


2 

















2 


1 


2 


1 


2 





1 





2 





2 




1 





1 


2 


2 


2 











2 











1 














2 


2 





2 





2 


1 








2 


2 








2 





2 





1 


2 








1 


1 








2 


2 




1 


2 








2 







2 








2 


1 


V o 


1 


1 


1 


1 


1 


1 


1 


1 


1 


o / 



In particular, note that each row and column sums to 9 while each broken diagonal 
sums to 8. 



4.4. Qualitative behavior of eigenvalues. By the triangle inequality one ob- 
tains the trivial bound \K(a, b)\ < p — 1 for all a, b. However, a significant amount 
of cancellation can occur in the sum ( |1.1[ ). The famous Weil bound asserts that 

\K(a,b)\<2 y /p, (4.15) 

whenever p \ ab [53]. A complete proof, based on Stepanov's method [21], can be 
found in the recent text [9l Thm. 11.11]. 

For a n x n real symmetric matrix X we let 

AoPQ < AipO < ■•• < Xn-l(X) 

denote the eigenvalues of X, repeated according to multiplicity. We are concerned 
here with the qualitative behavior of the eigenvalues of the (p + 2) x (p + 2) matrix 
T := Ti (3.12) and its p x p upper-left principal submatrix A := A\. 
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—2 /» 




"U v / 


A„ o(A) 


A„^i (T) 


2 fr> 


7 


—5.2915 


— Z.092UZ 


— z.oooy4 


Q OTQ1 1 
O.BYiill 


4.49390 


5.2915 


11 


-6.63325 


-5.71695 


5.3493 


4.48588 


4.79575 


6.63325 


29 


-10.7703 


-9.50028 


-9.43532 


8.89626 


9.06824 


10.7703 


71 


-16.8523 


-15.8699 


-15.8149 


14.1059 


14.1728 


16.8523 


113 


-21.2603 


-20.9713 


-20.8836 


19.5715 


19.6731 


21.2603 


229 


-30.2655 


-29.8296 


-29.75 


29.9351 


30.0001 


30.2655 


379 


-38.9358 


-38.2481 


-38.2008 


37.4232 


37.4756 


38.9358 


541 


-46.5188 


-46.4712 


-46.4221 


46.3519 


46.3885 


46.5188 


863 


-58.7537 


-58.5638 


-58.5258 


57.613 


57.6483 


58.7537 


1223 


-69.9428 


-67.6103 


-67.5843 


69.0147 


69.0451 


69.9428 


1583 


-79.5739 


-79.328 


79.3055 


77.3993 


77.4206 


79.5739 


1987 


-89.1516 


-88.7625 


-88.7417 


88.7745 


88.7849 


89.1516 



TABLE 4.3. The smallest and the second largest eigenvalues of the (p + 2) X 
(p + 2) matrix T and its (p — 1) X (p — 1) principal submatrix A. 



By the Weil bound (4.15) and a standard result relating the eigenvalues of a 



hermitian matrix to those of a principal submatrix [5[ Thm. 4.3.15] we have 



for < j < p — 2. In particular, it follows from (4.15) and the preceding chain of 
inequalities that 

-2y/p < M T ) < M A ) < A 2 (T) (4.16) 

and 

Ap-i(T) < A p _i(A) < A P+1 (T) < 2^P. (4.17) 

Using the Weil bound, we now write 

K(u) — 2^/pcos6 p (u) 

where p (u) E [0, 7r]. The vertical Sato-Tate law [TJ [11] states that as p — > oo 
the sequence of angles 9 p (u) becomes equidistributed with respect to the Sato- 
Tate measure p = ^ sin 2 9 dO on [0, 7r]. Thus for any fixed 5 > there are at 
least three values of 6 p (u) in each of the intervals [0,5] and [it — 5,w] when p is 
sufficiently large. In light of (4.16) and (4.17), we see that lim p _j. 00 Ao(A) 

2(A) 



and linip-j.oo X p 



2^p 
even 



2y/p. This behavior is clearly reflected in Table 4.3 
for relatively small values of p. 

The preceding argument relies upon a deep result of Katz [IT]. On the other 
hand, the matrix A is quite concrete and it enjoys many unusual combinatorial 
properties (Subsection 4.3). One might hope to estimate the eigenvalues of A 
directly to obtain an elementary proof of a Weil-type bound. The following result 
indicates that the error incurred using such an approach would not change the order 
of magnitude of the resulting estimate. 

Theorem 4.1. If p > 5 is an odd prime, then 

m&x{\K(u)\ : u = 0, 1, . . . ,p- 1} < max{|A (A)|, |A P _ 2 (A)|} + y/p^T. 

In other words, an estimate of the form max{|Ao(A)|, |A p _2(^4)|} < C^fp leads to a 
Weil-type estimate of the form tclwJ^K(u)\\ < C' ' s/p- 
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Proof. First note that if Aij are block matrices of the appropriate size, then 



An 



A 12 
A22 



A ln \ 
A 2 , 



< 



/Pull 
IIA21II 



U12W 
\\A22 1| 



\A ml \\ \\A m2 \\ 



\\Aln\\\ 

\\A 2n \\ 

\\AmnW J 



Moreover, if each is a nonnegative multiple of an all l's matrix, then equality 
holds. Therefore 



\\T-A®0, 



3x3| 



/- 







\ Vp~^ 











1 



1 



1 















000 

Vp^z 

y/p^l 000 



The theorem now follows from the triangle inequality for the operator norm. 
5. Applications to Kloosterman sums 



□ 



In the following, we employ the matrices T — T\ (3.12|, D — D\ (3.11), and U 



(3.131 constructed in Subsection 3.2 As before, we let A — A\ and B — B\ denote 



the upper-left p x p and (p — 1) x (j> — 1) principal submatrices of T. 



5.1. Basic Kloosterman identities. Using the character table of G (Table 2.2) 
we can derive a number of identities involving Kloosterman sums. For instance, 
taking the inner product of the first column with the (p + 2)th column yields 



p-i 

u=0 



(5.1) 



In particular, this gives another proof of (4.13) since trB = trT = trD = f — 1 + 

If c ^ 0, then taking the inner product of the first and the cth of the columns of 
the character table leads to 



p-i 

K{u)K(cu) = -p. 

Taking the inner product of the first column with itself we find that 

p-i 



(5.2) 



u=0 



# 









1 


si) 










(I 















y, z € Z/pZ 
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where Cg(-) denotes the centralizer of an element of G. This yields the well-known 
formula [HI eq. 3.7]: 



p-i 



(5.3) 



Since the matrices D and T (given by (3.11) and (3.12), respectively) are unitarily 



similar, it follows that the sums of the squares of their entries must be equal (i.e., 
tr£> 2 = trT 2 ). Thus 

(P 2 - p) + (-1) 2 + (P ~ I) 2 = tr B 2 + 2{p - 1) + 4(p - 2) + 2 (5.4) 



trD 2 



trT 2 



by (5.3). The preceding also yields tr_B = 2p — 9p + 10, which provides another 



proof of (4.11). In fact, since we already have an independent proof of (4.11), we 



could work backward from (5.4) to provide another proof of (5.3) 



From (5.2) and (5.3) we obtain (for c ^ 0, 1) 
p-i 



J2[K(u)-K(cu)] 2 = 2p\ 

ii=0 

p-i 

[K{u) + K{cu)} 2 = 2p 2 - 4p. 



(5.5) 
(5.6) 



Other such quadratic identities (e.g., [13 eqs. 3.5, 3.8]) might be deduced from 
Table 2.2 using the generalized orthogonality relations [3 Thm. 2.13]: 



Xi{GH) Xj (G-^) 



Gee 



' Xi(i) ■ 



Applying [10, Thm. 30.4] (see also [HI Prob. 3.9]) we obtain the formula 

2p 



IGI 



E 



X u (Gi)Xu(Gj)xu(Gk) 



|C G (G 4 )||C G (G,)|^ X u(I) 
which for i — 1 and 1 < j, k < p — 1 yields an identity equivalent to [13; Prop. 6.2]: 



p-i 



K{u)K{ju)K{ku) = [ /3( M fc ) ) / + . 2p 



u=0 



(this can also be easily deduced by computing the (J, k) entry of T = UDU). By 
( |4.2[ ) we have /3(1, 1, 1) = —3 from which we obtain [TU eq. 1], [15j eq. 3.22], and 
eq. (70)]: 



p-i 



p 2 + 2p ifp=l (mod 3), 
V " y [-p 2 + 2p Hp = 2 (mod 3). 



,*r d (u) = 



(5.7) 
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5.2. Quartic formulas and Kloosterman's bound. Computing the (j, j) entry 
of T 2 = UD 2 U, we obtain 



-J P "£k(u) 2 K(ju) 2 + 1+A 

P \u=0 ) 



entry of UD 2 U 



1 + 2(p-3) + f if j = 1, 

2 + 2(p-5)+2 if j^land(i) = l, 
2(p-3) + 2 if land (£) = -!, 



(j, j) entry of T obtained from Table |4.l| 

leading us to [T3l Prop. 6.3] and [HI eq. 3.18]: 



p-l 


V 


- 3p 2 - 3p 


if J = 1, 






3p 2 - 3p 


if j ^ 1 and (j) = 1, 


u=0 


P 3 - 


p 2 - 3p 


ifjV land (*) = -!. 



(5.8) 



Based upon this we obtain the following result of Kloosterman himself [12] 
Theorem 5.1 (Kloosterman). |if(u)| < 2ipi for all u. 



Proof. Let j = 1 in (5.8 1, observe that |if(u)| < 2p 3 , then take fourth roots. □ 



We remark that the simple proof above achieves a better constant (namely 2* 
in place of 3i - see also [HJ eq. 72]) than the recent proof in [S]. 

Although it is not clear whether one can obtain the Weil bound ( |4.15 1 using these 
methods, we have at least demonstrated that the unitary similarity A — UDU* 
encodes enough information about Kloosterman sums to obtain nontrivial results. 
Furthermore, we can establish that the exponent | appearing in the Weil bound 



cannot be improved. Indeed, from (5.3 1 and (5.8 1 we have 

p—i p—i 
2p 3 - 3p 2 - 3p = K ( u ) 4 ^ max{if(w) 2 } ^ K 2 (u) 



u=0 



u=0 



< ma,x{K(u) 2 }(p 2 - p) 



whence 



max{K(u) 2 } > 



2p 2 -3p-3 



p-l 

In other words, there exists some u such that 



= 2p - 2 



p — 5 
p~^l 



>2(p-l). 



\K{u)\ > V2(p-l)z. (5.9) 

5.3. Symmetric functions of Kloosterman sums. Recall that the coefficients 
Cj in the expansion 

det(X - XI) = c X n + ciA™" 1 + • • ■ + c n 

of the characteristic polynomial of a n x n matrix X are given by Bocher's recursion 

1 



co = 1, 



[cj-itrJf + c,-_ 2 trX 2 



c tr X j ] 



Applying this procedure to the diagonal matrix X — diag(_Koi Ki, ■ ■ ■ , K P -i) and 
using (5.1), (5.3), (5.7), and (5.8), we obtain Cq = 1, Ci = 0, 

c 2 = i(tr 2 A-trA 2 ) 
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c 3 = -g [(trX) 3 + 2trX 3 - 3(trX)(tr X 2 )] 



p + 2 



c 4 = i [(trX) 4 -6(trX) 2 (trX 2 ) + 3(trX 2 ) 2 +8(trX)(trX 3 ) -6trX 4 ] 
= |p(p-3)(p 2 -3p-2). 
We therefore obtain 




(5.10) 



which agrees with [151 P- 403]. The preceding now yields formulas for certain 
symmetric functions of Kloosterman sums: 



E K 3 K k = -\{p 2 -p), 



0<j<k<p-l 



p + 2 



0<i<j<fe<p-l 

E KiKjK k Ki = |p(p - 3) (p 2 - 3p - 2). 

0<i<j</c<i<p-l 



A different approach to (5.10) can be based upon the fact that the matrix 



X = 



( 01,1,1 Cl,l,2 
Cl,2,l £1,2,2 



Cl,l,/ 
c l,2,/ 



c l,/,/ 



V 





1 1 



1 

1 











o j 



is similar to the diagonal matrix D = diag(iTi, K%, ... , Kf, — 1, —1. /) ( |3.11 1. In- 
deed, the first matrix is similar to the truncated matrix T (3.12), which is itself 



unitarily similar to D. Using the fact that one may add a multiple of one row 
(resp. column) to another inside a determinant, one easily obtains 13 Lem. 14]: 



Theorem 5.2. The Kloosterman sums K , K\, K 2l 
values of the matrix 



Kf are precisely the eigen- 



( Cl,l,l — / Cl,l,2 - / 
c l,2,l c l,2,2 



V" 



Cl,/,1 



Cl,/,2 



c l,l,/ - / 
C l,2,/ 

Cl,/,/ 







where the coefficients Cj are defined by (4.1 ) and f = p — 1. 



Before proceeding, we remark that modifications of our main construction apply 
to various generalizations of classical Kloosterman sums. For instance, one might 
consider the Galois field F p n in place of Z/pZ. Moreover, our general scheme also 
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applies to hyper-Kloosterman sums (the appropriate analogue of our group G is 
discussed in [T31 p. 16]). 

6. Ramanujan multigraphs 

Certain principal submatrices of the Tj can be used to construct multigraphs 
having desirable spectral properties. To be more specific, a multigraph is a graph 
that is permitted to have multiple edges and loopsQ 

Associated to a multigraph Q is its adjacency matrix A(Q), the real symmetric 
matrix whose rows and columns are indexed by the vertices v\ , v% , . . . , v „ of Q and 
whose (J, k) entry a^fc is the number of edges connecting Vk to Vj. In particular, if 
j = k then ajj counts the number of loops attached to the vertex Vj . We refer to 
the eigenvalues of A{Q) as the eigenvalues ofQ. 

The degree of a vertex is the number of edges terminating at that vertex. We 
say that a multigraph Q is d-regular if each vertex has degree d. In this case d is an 
eigenvalue of Q with corresponding eigenvector (1,1,..., 1). On the other hand, — d 
is an eigenvalue of Q if and only if Q is bipartite, in which case the multiplicity of — d 
corresponding to the number of connected components of Q. An easy application 
of the Gerschgorin disk theorem indicates that every eigenvalue of Q belongs to the 
interval [— d, d]. We therefore label the eigenvalues of a d-regular multigraph Q, 
according to their multiplicity, as follows: 

d > A o (0) > Xi(G) > ■ ■ ■ > \ n -i(G) > -d. 

The eigenvalues of Q of a (i-regular multigraph which lie in the open interval (— d, d) 
are called the nontrivial eigenvalues of Q . We let X(G) denote the absolute value of 
the nontrivial eigenvalue of Q which is largest in magnitude. 

Following [20] , we say that a Ramanujan multigraph is a d-regular multigraph Q 
satisfying 

\{Q) < 2Vd-l. 

For instance, the Petersen graph is an example of a 3-regular Ramanujan graph 
(see Figure [T]). There is a vast literature dedicated to the study of simple (i.e., no 



FIGURE 1 . The Petersen graph is 3-regular and has characteristic polyno- 
mial (z — 3)(z + 2) 4 (z — l) 5 . The nontrivial eigenvalues 1 and —2 are both 
smaller than 2\/2 in absolute value whence the Petersen graph is Ramanujan. 

loops or multiple edges) Ramanujan graphs. We refer the reader to the seminal 
papers [THl HH1 US] and the texts [3] and [T7] for more information. 



The terminology in the literature is somewhat inconsistent. The term multigraph is some- 
times reserved for graphs with multiple edges but no loops. If loops arc present, then the term 
pseudograph is used. 
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We are now in a position to construct a family of Ramanujan multigraphs: 



Theorem 6.1. Let p > 5 be an odd prime, 1 < i < p — 1, and let (3(i, j, k) be given 
by $L2§. 

(1) The multigraph Q whose adjacency matrix is given by the matrix 



a j,k = 1 

is a (p — 2) -regular Ramanujan multigraph on p vertices. 
(2) If p = 3 (mod 4), then setting a^i = 1 and 

a j,k = 1 

otherwise yields a (p — 3) -regular Ramanujan multigraph on p — 1 vertices. 



Proof. The regularity of the resulting multigraphs are ensured by (4.8 1 and the 
general form (3.121 of T^. The fact that they are Ramanujan follows from (4.16) 
and ( |4.17| . □ 

For instance, letting p = 7 and i = 1 we obtain the adjacency matrix 

/202 1000\ 



10 10 2 1 
2 2 1 

1 1 2 1 
2 2 1 
2 2 1 

\ 1 1 1 1 1 / 



(6.1) 



corresponding to the multigraph depicted in Figure [2] 



FIGURE 2. The Ramanujan multigraph corresponding to the adjacency ma- 
trix \6.1\ . 
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